arXiv:1507.04660v2 [math.PR] 21 Jan 2016 


The Vertex Reinforced Jump Process and a Random 
Schrodinger operator on finite graphs 
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Abstract 

We introduce a new exponential family of probability distributions, which can be viewed as 
a multivariate generalization of the Inverse Gaussian distribution. Considered as the potential 
of a random Schrodinger operator, this exponential family is related to the random field that 
gives the mixing measure of the Vertex Reinforced Jump Process (VRJP), and hence to the 
mixing measure of the Edge Reinforced Random Walk (ERRW), the so-called magic formula. In 
particular, it yields by direct computation the value of the normalizing constants of these mixing 
measures, which solves a question raised by Diaconis. The results of this paper are instrumental 
in[T6l, where several properties of the VRJP and the ERRW are proved, in particular a functional 
central limit theorem in transient regimes, and recurrence of the 2-dimensional ERRW. 


1 Introduction 

In this paper we introduce a new multivariate exponential family, which is a multivariate generalization 
of the inverse Gaussian law. This exponential family is associated to a network of conductances and 
provides a random field on the vertices of the network, the latter having the remarkable property that 
the marginals have inverse gaussian law and that the field is decorrelated at distance two. 

This exponential family is mainly motivated by the study of two self-interacting processes, namely 
the Edge Reinforced Random Walk (ERRW) and the closely related Vertex Reinforced Jump Process 
(VRJP), but it could also find some applications in other topics, such as Bayesian statistics for 
instance. Note that Diaconis and Rolles introduced in 2006 a family of Bayesian priors for 
reversible Markov chains, similarly associated to the limit measure of the ERRW. 

More precisely, we consider a non-directed finite graph Q = {V, E) with strictly positive con¬ 
ductances Wij = Wj^i on the edges. Denote by the discrete Laplace operator associated to 
the conductance network (Wjj) and write Wi = J2j{i exponential family provides a 
random vector of positive reals {f3j)j^v such that 

Hp := -A^ + 

is a.s. a positive operator, where V = 2/3 — W is the operator of multiplication by {2(3i — Wi) and 
2/3 — W is considered as a random potential. We prove in Theorem [3] that if the Green function is 
defined by G = (W^)“\ then the field (e^-’) giving the mixing measure of the VRJP starting from 
to, c. f. m, is equal in law to (G(fo, j)/G(foGo))- 

This has several consequences. Firstly, it relates the VRJP to a random Schrodinger operator 
with an explicit random potential with decorrelation at distance 2. Note that Anderson localization 
was the main motivation in the papers of Disertori, Spencer, Zirnbauer ([21 HH]): in these works the 
supersymmetric field related to the mixing measure of the VRJP (c.f. [T4|) is viewed as a toy model 
for some supersymmetric fields that appears in the physics literature in connection with random band 
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matrices. Secondly, it enables one to couple the mixing fields of the VRJP starting from different 
points. Finally, using the link between VRJP and ERRW [H], it yields an answer to an old question 
of Diaconis about the direct computation of the normalizing constant of the ‘magic formula’ for the 
mixing measure of ERRW. 

Results of this paper are instrumental in [H], where the representation in terms of a random 
Schrodinger operator is extended to infinite graphs. Interesting new phenomena appear in the tran¬ 
sient case, where a generalized eigenfunction of the Schrodinger operator is involved in the represen¬ 
tation. Several consequences follow on the behavior of the VRJP and the ERRW in jl5]: in particular 
a functional central limit theorem is proved for the VRJP and the ERRW in dimension d >3 at weak 
reinforcement, and recurrence of the 2-dimensional ERRW is shown, giving a full answer to an old 
question of Diaconis. 

The paper is organized as follows. In Section [2 we define the new exponential family of distribu¬ 
tions and give its first properties. In section [2 we discuss the link between the exponential family and 
the Vertex reinforced jump processes. In Section |4]we consider the ERRW and answer the question 
of Diaconis. Sections |5] and 12 provide the proof of the two main results, namely Theorem [1] and 
Theorem |2 


2 A new exponential family 

Let V = {1,..., n} be a finite set, and let be a set of non-negative reals with Wij = 

Wj^i > 0. Denote by E the edges associated to the positive Wij, i.e. consider the graph Q = {V, E) 
with {i,j} G i? if and only if Wij > 0, and write i ~ j if {i,j} E E. Let Ag be the graph distance 
on Q. 

When A is a symmetric operator on (also be considered as a V xV matrix), write A > 0 if 
A is positive definite, and |y4| for its determinant. 

Theorem 1. Let P = {Pi,j)i<i,j<n be the symmetric matrix given by 


P = 


0 i = j, 

W,, i^j. 


For any 0 E M”, we have 




dfd 


exp ( - I ■ n ^ 

\ H.iteE / i=l V * 


where df3 = dfdi ■ ■ ■ djdn, snd 213 — P is the operator on defined by 

1(2/3 - P)f]{i) = 2/3,/(*) - WtjfU). 


( 1 ) 




Definition 1. The exponential family of random probability measures v^’^{dj3) is defined by 


,Wfii 


m = ^20-p>o(^T'^^w\-{o,i3)+ Y 


{hj}CE 


rij 


:d(3 


where {9, f3) = simply write for in the case where 9i = 1 for all i eV. 
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The proof of Theorem [T] is given in Section [5l We deduce the following simple but important 
properties of the measure . 

Proposition 1. The Laplace transform is 


= exp 


(^v^i+0iV^j+~ I ■ n 

{i,j}eE ) i=i 




e, 

Aj + 9i 


Moreover, if (3 is a random vector with distribution , then 


The marginals fdi 

^ =, 1 ) 


are such that 




is an Inverse Gaussian distribution with parameters 


• lfVi<zV,V 2 <zV are two subsets ofV such that dgiVi, V 2 ) > 2, then {I3i)i^vi {{^j)jcV 2 
are independent. 


Proof The Laplace transform of can be computed directly from Theorem [H from which we 
deduce independence at distance at least 2. We can also deduce, by identification of the Laplace 
transforms, that the marginals of this law are reciprocal inverse gaussian up to a multiplicative 
constant. □ 


The family can be reduced to the case 6 = 1 by changing W , as shown in the next corollary. 

Corollary 1. Let {(3j)j^v be distributed according to Then {9f3) is distributed according to 

where 

It is clear from the expression of the Laplace transform that if the graph has several connected 
components then the random field ((3j)j^v splits accordingly into independent random subvectors. 
Therefore, we will always assume in the sequel that the graph G is connected. 


3 Link with the Vertex reinforced Jump process 


3.1 Vertex Reinforced Jump Process: definition and main properties 

In this section we explain the link between the exponential family of Section [2] and the Vertex 
reinforced Jump Process (VRJP), which is a linearly reinforced process in continuous time, defined 
in f5], investigated on trees in |3], and on general graphs by the first two authors in [14]. Consider 
as in the previous section a conductance network (Wij) and the associated graph G = {V,E). Fix 
also some positive parameters (0j)jgv' on the vertices. Assume that the graph G is connected. 

We call VRJP with conductances (Wij) and initial local time {(pi) the continuous-time process 
{Xt)t>o on V, starting at time 0 at some vertex Iq eV and such that, if F is at a vertex i G V at 
time t, then, conditionally on {Ys, s <t), the process jumps to a neighbour j of i at rate 


where 


Lj{t) := (pj + / l{Ys=j}ds. 


3 










The following time change, introduced in El. plays a central role. Let 


D(t) = (L?(«) - 4^), 

i&V 


( 2 ) 


define Zt as the time changed process 

Zt = YD-^{t)- 

Let be the local time of Z at time t (that is, ij{t) = lz^=jds). Conditionally on the past, 

at time t, the process Z jumps from Zt = i to a neighbour j at rate (c.f. ng, Lemma 3) 


2 


We state below one of the main results of [T4j, Proposition 1 and Theorem 2. The theorem was 
stated in ng in the case 0 = 1, this version of the theorem can be deduced by a simple change of 
time, details are given in Appendix [Bl 


Theorem 2. Assume that Q is finite. Suppose that the VRJP starts at iq. The limit 

exists a.s. and, conditionally on U, Z is a Markov jump processes with jump rate from i to j 




Moreover (Uj) has the following distribution on {(««), = 0} 

n 


Q^’^idu) = e~ du, (3) 


with du = njGy\{io} 


D{W,u) = J2 n 

T {i,j}GT 

where the sum runs on the set of spanning trees T ofQ. We simply write for 

The fact that the total mass of the measure is 1 is both a non-trivial and a useful fact: in 

*0 _ _ _ 

particular, it plays a central role in the delocalization and localization results of [3 [in]. In [H] it 
is a consequence of the fact that it is the probability distribution of the limit random variables U. 
In [lU] it is proved using a sophisticated supersymmetric argument, the so-called localization principle. 
Theorem [3] below provides a direct ’computational’ proof of that result, based on the identity ([1]) and 
on the change of variable in Proposition [2] that relates the field (uj) to the random vector (0^) in 
Definition [H 
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3.2 Link with the random potential f5 

The second main result of this paper enables us to construct the mixing field e“ defined in the previous 
subsection from the random potential {f3j) defined in Definition [H It gives also a natural way to 
couple the mixing measure of VRJP starting from different points. 

Let us first state the following Proposition [21 which provides some elementary observations on 
the Green function. 

Define 

V = {{(], 2 / 3 -P> 0 }. 

Proposition 2. Let jS & V, and let G be the inverse of {2(3 — P). Then {G{i,j)) has positive 
coefficients. Define {u{i, j))i,jcv by 

G{z,z)- 

Then for G V, the function j —)■ u{io,j) is the unique solution j i—)■ Uj of the equation 

f = (3i, io 

[Mio = 0, 


In particular {u{io, j))jcv is {f^j)jcv\{io} measurable. Moreover, at the site io we have 


Ao 


1 

2 G{io, io) 


+ E 


Theorem 3. Let (3 be a random potential with distribution iy^’f'^{d(3) as in Definition \I\ and let 
{u{i, j))i,j€V be defined as in Proposition\2i Then the following properties hold: 

i) The random field {u{io, j))jev has the distribution of the mixing measure Q^’‘^{du) of the VRJP 
starting from io with initial local time {((iji^y 

a) The random variable G{io, io) has the distribution of 1/(27), where 7 is a gamma random variable 
with parameters (1/2,1/0^^). Moreover, G{io,io) is independent of {(3 j)j^ig, and thus also of 
the field {u{io,j))jev- 

The proofs of Proposition [2] and Theorem [3]are given in Sectional The next Corollary [2] describes 
how to construct the random potential (3 from the field u of Theorem [21 

Corollary 2. Consider a VRJP with edge weight (Wij) and initial local time ((/)j)jg\/- starting at io. 
Let {ui)i^v be distributed according to of Theorem [21 Let 

ft = ^ E (5) 

j- 3 ~i 

Let 'y be a Gamma distributed random variable with parameters (i, 1/0^^), independent of {uj), 
and let 

/3i = Pi + lio7. (6) 

Then (3 has the law of Definition [II 
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Corollary [2] indeed follows directly from Theorem [3] and Proposition [21 the law of in (|6]) is 
uniquely determined by the laws of snd 7 independent from hence it is sufficient to 

show that, if /9 has distribution {d(3) and u is defined from (|4]) by Proposition [2 then {ui)i^v 
indeed has distribution , and 7 = — Ag = l/(2G(fo, * 0 )) has distribution r(l/2, 

which follows from Theorem 12 

As mentioned in the introduction, Theorem [2 has several consequences. Firstly it explicitly 
relates the VRJP to the random Schrodinger operator —+ V, where V is the random potential 
Vi = 2{3i — Wi. Secondly it yields a natural coupling between the random fields {uj)j^v associated 
with the VRJP starting from different sites, since the exponential family {fdiji^v gives the same role 
to each vertex of the graph, and {u{i, j))ij^v arises from these random variables {I3i)i^v Finally 
it also gives a computational proof of the identity j Q^'^{du) = 1, for any 9, as a consequence of 
Theorem H] that allows to define (dP) as a probability measure. 


4 Link with the Edge reinforced random walk and a question 
of Diaconis 


4.1 Definition and magic formula 

The Edge Reinforced Random Walk (ERRW) is a famous discrete time process introduced in 1986 
by Coppersmith and Diaconis, [4]. 

Let {aij){ij}^E be a set of positive weights on the edges of the graph Q. Let be a 

random process that takes values in V, and let J^n = cr{Xo,... ,X„) be the filtration of its past. For 
any e G i?, n G N, let 

n 

■^n(c) = fle + 'y ^ l{{Xfc_i,Xi;}=e} (7) 

k=l 

be the number of crosses of the (non-directed) edge e up to time n plus the initial weight Og. 

Then (Xn)nm is called Edge Reinforced Random Walk (ERRW) with starting point fo G V and 
weights {ae)eeE, if = *0 snd, for all n G N, 


IP(^n+l — j \ n) 




Zn{{Xn,3}) 
Zn{{Xn-, k}) 


( 8 ) 


We denote by law of the ERRW starting from the initial vertex and initial weights 

(a). 

A fundamental property of the ERRW, stated in the next theorem, is that on finite graphs the 
ERRW is a mixture of reversible Markov chains, and the mixing measure can be determined explicitly 
(the so-called Coppersmith-Diaconis measure, or ‘magic formula'). It is a consequence of a de Finetti 
theorem for Markov chains due to Diaconis and Freedman |7], and the explicit determination of the law 
is due to Diaconis and Coppersmith [41ITT11I2]- It has also applications in Bayesian statistics [2[IllH]- 

Theorem 4. flJ]/ 

Assume that Q = (V, E) is a finite graph and set a* = J2j {i i E V. Fix an edge 

Co incident to and define FLeo = {y ■- E E, > 0, ye^ = 1} (similarly let yi = Y^iaeVe)- 
Consider the following positive measure defined on FLeo defined by its density 

Mtjidy) ^ vw) n —■ ^ 

it 
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with 


T eGT 

where the sum runs on the set of spanning trees T of Q, and with 


C{a,io) 


n.gy r(i(ai + 1 - 


|y|-i 


rieGE ^(®e) 


Then is a probability measure on TLeo, ^nd it is the mixing measure of the ERRWstarting from 
io, more precisely 

J.ERRW,M^,) _ r pif(.)dM\:\y), 

where denotes the reversible Markov chain starting at io with conductances (y). 


4.2 The question of Diaconis 

The fact that Ai\'^\dy) is a probability measure is a consequence of the fact that it is the mixing 
measure of the ERRW. In fact it is obtained as the limit distribution of the normalized occupation 
time of the edges m 

f Zn{e) \ law ^{a) 

\Zn{eo)JaeE 

One question raised by Diaconis is the following 

(Q) Prove by direct computation that J Ai\'^\dy) = 1. (10) 

An answer was proposed by Diaconis and Stong |Ej in the case of the triangle, using a subtle change 
of variables. Also note that MerkI and Rolles offered in [12] analytic tools for the computation of 
the ratio of the normalizing constants of the magic formula for two initial weights differing by integer 
values, which may possibly be extended to provide the normalizing constant. 

We provide below an answer to that question. A first simplification comes from [T4], where an 
explicit link was made between the VRJP and the ERRW. 

Theorem 5 (Theorem 1, [H]). Consider (y„) the discrete time process associated with the VRJP 
(Yt) (i.e. taken at jump times) with conductances (Wij) and cj) = 1. Take now the conductances 
{WejeeE 3S independent random variables with gamma distribution with parameters {aejecE- Then 
the 'annealed' law ofYn (i.e. the law after taking expectation with respect to the random {We}) is 
the law of the ERRW (X„) with initial weights {ae)ecE. 

This immediately implies an identity between the mixing measures and indeed, by 

Theorem [21 (Yn) is a mixture of Markov jump processes with conductances which implies 

that for all 0-homogeneous bounded test functions (j) (i.e. (j){\y) = (j>{y), VA > 0), we have 

/ •P{i.Ve))Mt^[dy)= ( n "'* ( / dW. (11) 

JHeo eCE ^ ^ 

with dW = riees'^^e- This identity was checked by direct computation in section 5 of [14] . 
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Now, the fact that f Qj^(du) = 1 is a consequence of the computation of the integral ([I]) i 


in 


Theorem [Hand the change of variables in Theorem O as explained at the end of Section [31 Therefore 


/ dMl(y) = 1. 

dyeQ=l 

Note that this fact can be used to prove directly that M.^^[dy) is the mixing measure of the 
ERRW starting from initial condition (a) and initial vertex iq. Indeed, for any finite path a : fo 

let N{i) (resp. N{e)) be the number of times vertex i (resp. edge e) is visited (resp. 

crossed): 

N{i) = \{k] 0 < k < n — 1, ik = i}\ 

N{e) = \{k; 0 <k <n-l, {4,4+i} = e}|. 

The probability of a for the reversible Markov chain of conductance y is 

n 

lleS-B ye 




riiGy Vi 


N(i) 


The integration ofp^^((T) w.r.t. (?/) can be computed by changing the constant r(ae) to r(ae + 
Ne) and r(i(aj + 1)) to r(|(ai + 1) + Ni). Using the property r(x + 1) = a:r(x) and the notation 
(a, n) = nfc=o(® + ^)’ deduce 


p\M)dMUy) = 


which is the probability of an ERRW to follow the path a. 


5 Proof of Theorem U 

Lemma 1. Let P = iPi,j)i<i,j<n be a symmetric matrix with 


P = 


0, i = j, 


and let (d be a diagonal matrix with entries fdi, i = 1,... ,n, such that M = 2(d — P is positive 
definite. 

Let L be the lower triangular n x n matrix and U be the upper unitary (with 1 on the diagonal) 
upper triangular matrix such that M = LU (i.e. the LU decomposition of M), which exist and are 
unique. 

Then 


U = 


where {xi}i<i<n 3nd (Wjj)i<j<j<„ are defined recursively by 


/ Xi 



\ 

0 

X2 

-H2,n 




Hn—l,r 

i 

VO 

0 

Xn 

/ 


'Hia = ^1,. 

= W".. + Er=\ 




^ 1 

^k 


j > 1 

i > 2, j >i 
i > 1. 
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Furthermore, 


Xi = 


M{1,... .. ,i) 


where M{I\J) is the minor of matrix M that corresponds to the rows with index in I and columns 
with index in J. 

The result follows directly from (2.6) of [IT], but we prove it in Appendix O for completeness’ 
sake. 

Claim 1. For any 6i > 0, 6 * 2 > 0 , 

Proof The case 6*2 = 0 corresponds to the normalisation of the r(i) variable. The case 6*2 > 0 
corresponds to the normalization of the Inverse Gaussian law IG(|^, 


^2 ’ ^2 ^ ' 


□ 


Let us now prove Theorem [H In the sequel we take the convention, given any real sequence 

(afc)feGN. that flfc = 0 if n > m. 

By Lemma [H 


n k—1 tj2 

,Xk 




k=l 


2 2xi 

k=l 1=1 ‘ 


i=l 


^ + J_( V e,Hf 

2xj ^ 


k=l+l 


Define 


^ ; 


.+ \ { 0 })" 
iXi)l<i<n 


V 


i-l tt2 
Xi 


k=l 


l<i<n 


Then \l/ is a bijection, since a symmetric matrix is positive definite if and only if all of its diagonal 
minors are positive. Its Jacobian is 2“”, hence it is a diffeomorphism. 

Therefore 


/:= / 1 


{2,a-p>o}' 


exp(-0/3) 


Let, for all 1 < / < m < n, 


df3 = 


/ -(-2 


k=l+l 


y/Xi ■■■Xn2^‘ 


-dx. 


E li'j'/si] + E 

\j = 771+1 

Rl,7n 


k=l+l 


Sim — 


+ 


2xi 


Note that Ri^m (resp. Si^m.) only depends on xi, ... xi-i (resp. xi, ... xi). 
Let, for all 1 < m < n. 


dm. • 


Tr- V 1=1 J 


dxi ■ ■ ■ dXr, 


■■■Xr. 


We will take the convention that, if m = 0, the integral of dxi ■ ■ ■ dxm is 1 , so that Jq = 1. 
Note that / = In/2'^. We also have the following lemma. 
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Lemma 2. For all 1 < m < n, we have 


■^m xl Q 6Xp j ^ ^ \/| ^m—1- 

™ ' j=m+l 


Proof. Using Claim [H we deduce 

Im= exp 


m —1 


Rm,m \ q 

+ ^ 

^Jvr, 


Lm 


1=1 


dxi ■ ■ ■ dXr. 

= f exp (-./itza; - E s,„) 

Jk--! \ ^ J • • • Xm-1 

Now R,n,m = (EJ=m+i and 

m—1 

Hm,j = 




l=l 


Xl 


so that 


n m—1 j— — n 

x/Rm^mOm = ^ 

l'=m+l i=l ?'=m+l 


j=rn+l 

On the other hand, for all 1 < / < m — 1 , 


C C _ JJ 

Ol,m — Ol,m-l — -- 2_^ 


Xl 


ji =772+1 


m—1 


Therefore 

+ Rm,mdm + Si^m = ^m,j + * 5 '/ 

Z=1 j=m+l 

which enables to conclude by (ED- 

We deduce from Lemma [21 by induction, that 


m—1 


1=1 


/,m —1 ? 


/=^ = lJJM^exp(- 

UjyeE 


2 ^ 2^ V + • • • 

which enables us to conclude. 


( 12 ) 


□ 


6 Proof of Proposition |2l and Theorem 31 

6.1 Proof of Proposition [2] 

Fix io E V, and let (3 E V. Let us first justify the existence and uniqueness of M(fo,f) defined by 
the linear system (01). As (2/3 — P) is an M-matrix, its inverse G satisfies G{i,j) > 0 for any i,j- A 
solution (uj) of equation ( 0 ]) is necessarily of the form = 2'yG{io,j) foi' some constant 7 G M. 
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The normalization Ui^ = 0 implies 7 = 2 G(J^~i^- Hence the unique solution of the system (S]) is given 
by Uj = u{io,j) defined in Theorem [3l 
Consider the following map: 

: V ^ {{uj)j^v ^ = 0} X (M+ \ {0}) 

(13) ^((mj),7 ), (13) 

where (uj) is the unique solution of the system ([4]) and 7 = 

We first prove that <l> is a diffeomorphism. By the previous argument it is well-defined and 
injective. Reciprocally, starting from ({uj),'y) on the right hand side, we define (/3j) by 

ft = E (14) 

It is clear that with this definition, (uj) is the solution of (|4]) with {(3j). It remains to prove that 
2(3 — P > 0\ it is a consequence Theorem (2.3)- (J30) of [T3] : 

Proposition 3. Let A & Zn = {M G M„(M), < 0, ifi 7 ^ j}. Then A is positive stabl^ if and 

only if there exists ■C > d§ with > cH and 

k 

^ ayjij >0, /c = 1,..., n. (15) 

We will choose a bijection a between V and {1,..., |f/|}, and apply Proposition Owith 

A = ((2/3 - -P)(T~i(i),o-io))i<j,i<|y|, ^ = (e“'"“^(ft)i<i<|\/|. 

Obviously, ^ S> 0, and > 0 follows from (2/3 —P)e“' = Sig/(2G(io,io))- Now fix any spanning 
tree T of the graph and its corresponding distance d on V throughout the tree. Choose a so that 
cr(fo) = \y\, and a(i) < a(j) if d(io,i) > d(io,j)- this implies that, for all k < \V\, there exists 
I > k such that fT/r-i(fe),< 7 - 1(0 > 0 and therefore that (flSi) holds. We conclude that 2/3 — P > 0. 

6.2 Proof of Theorem S] 

We give two proofs. 

First proof: We make the change of variable given by in ([T3|) and we now prove that if /3 
has distribution , then (u, 7 ) = <l)“^(/3) has distribution (8) r(i, ^). 

0 ^iQ 

Let J be the Jacobian matrix of (i.e. Jtj = J 7 ^ io Ji,io = ^), then 

{ ^i,io j ^0) 

-(3i \f i = j ^ io- 

We can factorize the fth row of J by for each i, then expand the resulting matrix according to 
the ioth column, and we find that 


^All of its eigenvalues have positive real part. 
^ rj means for any coordinate i, Li > 

> 0 means > 0 and ^ yf 0 
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On the other hand, by (fl^ we deduce 

|2/3-P| = 2'ye-^^^^^D{W,u). 


Let '0 be a positive test function. We have 






—S^ 7 
e ^0 ' 




= / <fe“,7)SPW — 


—cb^ 7 

e *0' 


7r7 


(^7. 


This concludes the proof of Theorem [3] and of Corollary [21 

Second proof: This proof does not make use of the explicit expression of law of U in ([3]), 
but rather deduces its Laplace transfom from direct computation of the probability of a path. Note 
that compared to the first proof, this one uses the representation of the VRJP as a mixture of Markov 
Jump Processes, cf Theorem 2 of [14] or Theorem [3] in section [31 and hence it uses implicitly that 
the measure is a probability measure. 

We will show that, if (m, 7 ) has distribution ( 8 ) r(i, ^), then fd = <I>(m, 7 ) has distribution 
^ which clearly implies the result. 

It follows by direct computation (see da. proof of Theorem 3) that the probability that, at time 
t, the VRJP Z has followed a path Zq = xq, Xi . Zt = Xn with jump times respectively in 

[ti, ti + dti], i = 1.. .n, where to = Q < ti < ... < tn < t = tn+i, is Ptdt, where 



1 

dt = —^xi-ixi dti, 
i=l 


with {ii)i^v = {^iit))iev local time at time t. 

On the other hand, using that, conditionally on P = {Uiji^v in Theorem [21 Z is a Markov Jump 
process with Jump rate /2 from i to j, this probability of a path is also qtdt, where 


Qt = 




and /3 is defined in ([5]). 

Let r = r(|,^). By identification ofp^ and qt we deduce that 



which shows that the distribution 0 r(i, has the same Laplace transform as in 

Proposition [H 
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A Proof of Lemma U 


Proof. We perform successive Gauss elimination on M to make it upper triangular. Denote by 
li,... ,ln the n rows of any n x n matrix. Firstly, let 


= m = 

( 

o-(l) 

-"1,2 

-"1,2 ■ ■ 

X2 

1 1 


\ ^l,n 

TTP) 

-"n,2 ■ ■ 

Xn'^ j 


where we set, for any 1 < i,j < n, = 2/3* and = Wij. 
We define a sequence of matrices recursively, such that 


/^(i) 


= 


^1,2 

X2 

0 


Hi'l \ 


-H. 


( 2 ) 

2,3 


-H. 


( 2 ) 

2,n 


X 


(fc-i) 

k-l 


-H. 


(k-i) 

k—l,k 


{k) 

^k,k+l 


-H, 


(k) 


k,k+l 


_ Tri^—P 
^k-l,n 

_ U'(^) 

^k,n 


V 0 


-w, 


(fc) 


k,n 


-H, 


(k) 


n—l,n 


^ n—l,n 


by the following rule: is constructed from by addition of columns Z^+i ^ 4+i + 


^^4, ...,4^4 + -^4 in In other words. 


where = 


H 


0 


1= J 
(k) 

i> j = k 

Otherwise 




It is easy to check that satisfy the following induction rule: 

) f, j > ^ + 1, 


r(fc)^ 


^(fc+1) _ TT{k) . t^k.i^k,- 


™(fc+l) _ Jk) _ 


i > k + 1. 


At step n, we have 

Tn-i ■ ■ ■ TiM = = 


(x? -h[] 
0 


( 1 ) 
2 

,^(2)' 

X2 


l,n 
(n) 

Xn 




\ 0 0 

Hence, it gives the LU-decomposition of M where L~^ = T = T„_iT „_2 • ■ - Ti and U = It is 

easy to check that 

Xi = x]' t = 1,... ,n 
i<j 

satisfy the recursion in the statement, and that x* = M(l,..., i|l, • • •, ..., / —1|1, • • • 4“!)■ 

□ 
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B Time rescaling 


Let Ys be the VRJP with conductances (W) and initial local time (0j)jgy defined in Section O Recall 
that Li{t) = (j)i + Jq Consider the increasing functional 74(s) = “ 1)’ 

time-changed process = ^ 4 - 1 ( 5 )- Let 


LAs) — 1 -l- 


^{Y=i}dS- 


We always denote by s the time scale of F, we can write 


s = A(s), ds 



1 


LAs). 


Obviously, F is a VRJP with edge weight WijAiAj ^nd initial local local time 1 : that is, conditionally 
on F jumps from i to j at rate 

^^ijAiAjLj (s). 

Note for simplicity 

WA = W,,ct>^Ay 

We can apply na Theorem 2 to F. Let 


D(S) = 5^L,(Sf-l, 

i 


and set Zi = F^-i(f), with local time £i(t) = proposition 1 of [T4] translated in time 

scale L (cf relation (2.1) of [H]), we have that logLi(s) — jfJ2jev^^sLj{s) converges a.s. when 
s —00 to a random vector with distribution given by (3.1) of theorem 1 of [14], where the weights 
{Wij) are replaced by (Wfj)- Changing to variables Ui ^ Ui — Ui^, we deduce 


lim logLi(s) - logZio(s) = Ui 

s—)-oo 


exists and has distribution 


«'(*) = 


v/2^ 


■N-l' 


=-E 


'j£V P 2 Yir. 


:W, 




and that Z is a mixture of Markov Jump Process with jumping rates We now come 

back to {Zt). Recall that Zt = where D{t) is defined in ((2]). ,!,From this we have 


i = D{A{D-\t))), 


and 


di 


1 Ly. (s) 
HLvsis) 



This implies that (Zt) is a mixture of Markov Jump processes with jumping rates log^i 

By simple change of variables, Ui J- log Ai — log Aio has distribution 




-T 




■N-l 


■jev “jg- 


-1 y. 
2 
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